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NGC 3603 is the most massive and luminous visible starburst region in the Galaxy. 
We present the first Chandra/ ACIS-I X-ray image and spectra of this dense, exotic 
object, accompanied by a deep cm-wavelength ATCA radio image at similar < 1" 
spatial resolution, and HST/ground-based optical data. At the S/N > 3 level, Chandra 
detects several hundred X-ray point sources (compared to the 3 distinct sources seen by 
ROSAT). At least 40 of these sources are definitely associated with optically identified 
cluster O and WR type members, but most are not. A diffuse X-ray component is also 
seen out to ~ 2 ; (4 pc) from the center, probably arising mainly from the large number of 
merging/colliding hot stellar winds and/or numerous faint cluster sources. The point- 
source X-ray fluxes generally increase with increasing bolometric brightnesses of the 
member O/WR stars, but with very large scatter. Some exceptionally bright stellar 
X-ray sources may be colliding wind binaries. The radio image shows (1) two resolved 
sources, one definitely non-thermal, in the cluster core near where the X-ray /optically 
brightest stars with the strongest stellar winds are located, (2) emission from all three 
known proplyd-like objects (with thermal and non-thermal components), and (3) many 
thermal sources in the peripheral regions of triggered star- formation. Overall, NGC 
3603 appears to be a somewhat younger and hotter, scaled-down version of typical 
starbursts found in other galaxies. 

Subject headings: open clusters and associations: individual (NGC 3603) — stars: early- 
type — stars: winds, outflows 


1. Introduction 

Hot stars, by means of the massive supersonic winds and the supernovae they produce, have 
a profound influence on their environment and on galactic and chemical evolution in the Universe. 
A sufficient density of massive stars can drive galactic-scale outflows, such as that seen in M82. 
While many advances have been made in the study of such starburst galaxies, they typically are 
too distant to study in great detail. Consequently, any local analogues, such as giant H II regions 
in our own Galaxy and in the Magellanic Clouds, are of great interest. 

Giant H II regions - particularly if relatively compact yet nearby - provide nearly co-eval sam- 
ples of large numbers of resolved massive stars. Deep X-ray imaging is an efficient way to measure 
the X-ray properties of a viable, homogeneous sample of OB/WR stars. We have therefore used the 
Chandra X-ray satellite to image NGC 3603 - the brightest, most luminous visible giant H II region 
with near solar abundance in the Galaxy, at an adopted distance of 7 kpc (Moffat 1983; Moffat et 
al. 1994; Drissen et all 995; Brandi et al.1999). The Chandra observations are complemented by 
cm-Avavelength radio imaging at the Australian Telescope Compact Array (ATCA) and by optical 
observations using the HST and groundbased telescopes. Such multi-wavelength data, all at similar 
spatial resolution, provide a clear advantage when it comes to understanding the physics of star- 
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burst regions. For example, X-ray and radio emission mechanisms are often related (e.g. in regions 
of shocked gas), and optical data are indispensable for identifying X-ray sources. This preliminary 
study will show how X-ray and radio properties of the cluster relate to the physical properties of 
the stars and gas in a starb irst region. 


2. Observations 

We have obtained a deep, 50-ksec, uninterrupted pointing at NGC 3603 with ACIST on Chan- 
dra in May 2000 (cycle AOl ), with a spatial resolution of ~ 1". The ACIS-I image was created from 
the standard photon events after they were cleaned to remove photons falling near node boundaries 
or near other bad pixel regions. Inspection of X-ray light-curves of a region centered on NGC 3603 
and on a background region to the north of the cluster indicate that there is no evidence for the 
presence of flares, so that no time intervals need to be excluded. The cleaned events were then 
spatially binned to produce an image of the central 3.6 / x at the full resolution of ACIS (0.492 
pixel" 1 ). We corrected the image for exposure variations across the detector and smoothed the 
resulting image with an adaptive filter. The smoothed, cleaned image is shown in Figure la. 

In addition to the X-ray data, we present radio data from four 12-hour shifts with the ATCA 
at 3 and 6 cm wavelengths taken in February, April, June and September 2000 (see Miicke et al. 
2001), resulting in a resolution of 1.0 ,/ x 0.8" at 3 cm with a noise level of ~ 0.15 mJy. The radio 
data use only the longest baselines (> 25 kA); this will tend to emphasize the small-scale structure. 
Finally we complement the X-ray and radio data by a broadband visual HST/WFPC image of 
the central core at 0.1" resolution, along with a 1" resolution image of the whole cluster and its 
surroundings (i.e. similar to the resolution of the X-ray and radio data) from the 0.9m telescope 
at CTIO. 

Fig. 1 shows overview linear-intensity scale images centered on NGC 3603, in (a) X-rays (0.5- 
10 keV, i.e. 24 - 1.2 A) and (b) groundbased medium-band optical light (A c = 0.465pm = 4650 A), 
both at the same size, scale and resolution (~ 1"). Also shown superposed (in blue contours) are 
the A = 3cm = 3 x 10 8 A radio data, which have been slightly adjusted in Right Ascension (RA) 
and Declination (DEC) so that the radio flux of the proplyd-like objects PI and P2 (P3 is outside 
the HST field in Ha) coincide best with their HST positions in Ha emission (see Miicke et al. 2001). 
Fig. 2 is an enlargement of the central field of Fig. 1, with superposed deconvolved ACIS-I X-ray 
(green) and smoothed HS I optical (red) images, and 3 cm radio contours. The Chandra image 
was shifted in RA and DEC by less than 0.5” to bring the brightest, most obvious HST sources in 
good positional agreement. Stellar coincidences between X-ray and optical appear in bright yellow. 
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3. Results and Discussion 
3.1. O-star X-ray Spectral Properties 

In principle it is possible to obtain X-ray spectra for bright detected sources directly from the 
ACIS data, since the CCD detectors measure the pulse-height amplitude for each detected photon 
from each source. In practice this is difficult to do directly since the brightest sources in the image 
suffer from pile-up (in which two or more photons arrive at a single pixel in one readout cycle and 
are counted as a single photon event), which causes a loss of flux and an increase in high-energy 
photons at the expense of low energy photons. Individual fainter sources which do not suffer from 
pile-up do not have enough counts for us to derive a well-defined X-ray spectrum. Assuming that 
stars of similar spectral type should possess similar X-ray spectra, we have defined a composite 
stellar X-ray spectrum by extracting spectra for individual faint sources of similar spectral types 
and co-adding the spectra to produce an unpiled spectrum of sufficient significance. There are 
seven member stars of early-0 spectral type in the core of NGC 3603 which are fairly well resolved 
in the ACIS image and detected at a signal-to-noise level >3, yet do not suffer from serious pile-up. 
These are H19, 49, 50, 52, 58, 61, 62 from Moffat et al (1994); see figure 1. 

Before extracting and coadding the spectra for these seven sources, we attempted to correct 
the detected X-ray photons for the effects of charge-transfer inefficiency (CTI), in which the pulse- 
height distribution changes as the charge migrates from the source position to the CCD readout. 
This is of concern here since the distortions due to CTI will be different for each of our sources 
depending on how far each source is from the chip readout. To minimize the effects of CTI on the 
extracted stellar X-ray spectrum we corrected the individual photon events using the algorithm of 
Townsley et al. (2000). 

We extracted the X-ray spectra from the cleaned, CTI corrected photons for these seven stars 
using extraction regions of about 1 // in radius. The composite spectrum is shown in Fig. 3. We 
extracted background from the cleaned, CTI-corrected photon event file in a nearby source-free 
circular region of radius 17". We binned the composite early-0 star spectrum so that each bin 
had a minimum of 10 counts and we then fit the binned, background-subtracted spectrum using a 
single-temperature Mewe, Kaastra, &; Liedahl (1995) emission model plus absorption by overlying 
cool material. The best fit for this model is shown superposed in Fig. 3. The derived temperature is 
kT = 0.85 ±0.1 keV and Nji = 1.3 ±0.1 x 10 22 cm -2 at the 90% confidence level. The fit using this 
single-component model is formally unacceptable (reduced x? ~ 1.29); most of the contribution to 
X 2 is produced by an apparent excess at E > 3 keV which is not matched by the single-component 
emission model. 12 The flux for this model is ~ 6 x 10 -14 ergs cm -2 s _1 in the energy range 0.5-10 
keV, suggesting that each star has an average X-ray luminosity of ~ 5 x 10 31 ergs s _1 for an assumed 
distance of 7 kpc. 


1 2 A multi-temperature model (e.g. as for three of the Orion Trapezium stars: Schulz et al. 2001) would undoubtedly 
fit better; however, we consider this to be over-analysis in the context of this overview paper. 



-5- 

3.2. Stellar Point-Source Members of NGC 3603 

The standard ACIS processing provides a source list produced by running a local cell detection 
algorithm (the celldetect tool from the Chandra Interactive Analysis of Observations [CIAO] version 
1.1) with the local signal- to noise ratio set to 3. This is probably adequate near the center of the 
field but may be less acceptable farther off-axis since sensitivity varies with off-axis position. The 
standard source list contains 174 sources. A casual inspection of the source positions compared to 
the adaptively smoothed exposure-corrected image suggests that the simple cell detection algorithm 
probably misses about 30-50% of the sources visible to the eye. We also analyzed the field using 
the wavelet source detection tool, wavedetect, available in CIAO version 2.2. The wavelet detection 
algorithm finds about 50 more sources compared to the local cell detect algorithm, although the 
X-ray fluxes from this are less reliable than with “celldetect”. As an alternative, we wrote our 
own simple cell-detection algorithm, which is similar, but not identical to the CIAO celldetect 
algorithm. First, our algorithm calculates the background level across the image iteratively omitting 
potential point sources. Then it assesses the significance of ~ point-like events using a specified 
statistical threshold (3<7 in our case), calculated after Gehrels (1986). Each point-like event is 
recognized as such after integrating the counts within a given distance (r = 2,3,4,.... pixels) from 
the potential source. Finally, the algorithm compiles a list of spatially-independent events (simply 
those which are separated by d > r). Using r — 4 pixels, we found 384 sources in the whole field. 
A more detailed analysis of the X-ray properties of detected X-ray sources, their identifications, 
and inclusion of upper limits for optical stars not detected in our Chandra pointing is currently 
underway; for the remainder of this paper we consider for analysis the source list derived from our 
own simple cell-detection algorithm, which yields the most sources. 

In Fig. 4 we show a plot; of X-ray (unabsorbed) versus bolometric luminosity of detected point 
sources labelled in Figs. 1 and 2. This sample contains all the sources for which X-ray fluxes 
are detected above the ~ 3-t level and for which optical magnitudes/colors/spectral-types reveal 
their membership as OB/WR stars in NGC 3603: cf. Melnick et al.(1989) and Skalkowski et al.(in 
prep.) for stars outside the dense core, and Moffat et al.(1994) and Drissen et al.(1995) for the 
cluster centre. Cluster membership is well established in the area of Fig. 1 (except possibly for 
some fainter stars in the core) down to V ~15 mag. This corresponds to L q 0 i ~ 8 x 10 38 erg s 1 , 
so that Fig. 4 is limited almost entirely by Lx- To derive Lsol f° r the member stars, we adopted 
the visual absolute magnitudes Mys from the above mentioned photometric studies. Then with 
constant bolometric correction {BC = My ~ Mq 0 i = 4.5 mag 13 ) and distance 7 kpc, we derived 
the M Bo is and Lb 0 /s. For conversion of ACIS count rate to flux L x we assumed a Raymond-Smith 
thermal spectrum with a temperature of 0.7 keV, typical of the X-ray temperatures of single massive 
stars (see, e.g. Corcoran et al. 1994) and consistent with other analyses (Chlebowski et al. 1989; 
Bergh“fer et al. 1997; this paper: see previous sub-section), and Nh derived using the relation 


l3 This is reasonable, since mosi. of the detected stars are of early O type or H-rich WNL and hence nearly constant 
BC 
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JVfi = 5.8 x 10 2l E(B — V ) (Bohlin, Savage and Drake 1978). Use of a Mewe, Kaastra, & Liedahl 
(1995) model instead of a Raymond-Smith model would lead to negligeable changes here. Note 
that an increase to kT =1.0 keV would increase the absorbed X-ray fluxes in the 0.5-10 keV band 
by ~15 %, which is relatively unimportant here. The unabsorbed (“dereddened”) X-ray fluxes are 
much more sensitive to the assumed spectral shape. However, further refinement is not justified 
with the present non piled-up spectra, which have only ~100 counts each at the best. We consider 
an apparent association between point X-ray and optical sources only if the separation does not 
exceed one resolution element, thus probably missing several potential associations (especially in 
Fig. 2). Neglecting this, we plot in Fig. 4 ~45 point source members of NGC 3603 with almost 
certain X-ray /optical associations. Of particular note are the following: 

(1) X-ray flux generally increases directly with bolometric flux in the usual way (e.g. Seward 
& Chlebowski 1982; Chlebowski et al.1989; Bergh "fer et al.1997), i.e. Lx ~ 10“ 7 LjB O /, extended 
here to higher bolometric luminosities. However, there is at least two orders of magnitude scatter 
in the X-ray flux around this simple scaling relation. The large scatter appears to be intrinsic to 
the NGC 3603 stars, since we only consider members of NGC 3603, with field stars excluded. 

(2) All of the brightest early-0 and all three of the WR stars are seen in X-rays. These 
relatively H-rich WR stars are by far the brightest optical sources in the cluster, all located in 
the core. They are probably main-sequence stars of extremely high luminosity with strong wind- 
produced emission lines, similar to the three bright H-rich WR stars in the R136 core of 30 Dor (De 
Koter et al.1997; Crowther & Dessart 1998). Thus, the 3 WR stars in NGC 3603 probably do not 
behave like genuine highly evolved (Galactic, single) WR stars of the WN sequence, for which Lx 
and LboI, ranging over 10 31 “ 34 and 10 38-39 - 5 erg s“ L , respectively, are uncorrelated (Wessolowski 
1996). 

(3) The two optically bright late- O/early- B supergiants Hl(OC9.7Ia) and H5(= Sher 25, 
B1.5Ia, with the SN1987A-like ring structures: see Brandner et al.1997), located ~ 20" outside 
the cluster core have no significant X-ray detection. This is likely a consequence of their relatively 
low Lb 0 i s combined with their slower winds. These two stars probably belong to a more extended, 
slightly older population, compared to the extreme youth of the rest of the known bright members 
of NGC 3603. 

(4) Several sources have relatively outstanding X-ray fluxes 14 for their LboIi in ~order of 
decreasing Lx/LboI • (a) The WR star C, which might be the suspected long-period (100s of days) 
WR binary of Moffat & Niemela (1984), possibly similar to the extremely X-ray bright H-rich WNL 
star WR25 in the Carina Nebula (Seward & Chlebowski 1982). (b) MTT 68 and MTT 71 (Melnick 
et al.1989; see Fig. 1), two heavily reddened ( Eb-v ~ 1-93, compared to the cluster core with 
Eb-v = 1.44) member OB stars located over 1' (2 pc projected separation) from the cluster core; 


14 Since no pile-up corrections were applied, the Lx -values for the X-ray brightest stars (C, 68, A1 + A2+A3) will 
be lower limits. However, even in the worst case, star C is estimated to be a factor ~ 2 too low in Lx, which is ~ 
negligeable in the context of Fig. 4 and its implications. 
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these two stars also stand out on i?OSAT X-ray PSPC images, (c) The WR star Al (although 
confused in X-rays with the 2 nearby but optically resolved 03 stars A2 and A3), which is now 
known to be a 3.771d binaiv (Moffat & Niemela 1984; Moffat et ah, in prep.), (d) The 03-4V 
+ 04-5V double-line system, star no. 11 in Fig.lb. Whether all or most of the enhanced X-ray 
emission arises from colliding winds in close binaries remains to be demonstrated. 

(5) Some remaining X-iay point sources have no obvious bright optical counterparts, yet their 
frequency increases towards the cluster center. Some of these might be correlated with optically 
faint but active pre-main sequence stars (e.g. Brandi et ah 1999). 

(6) Two significant but faint discrete radio sources above the noise level lie in the central 
cluster region of Fig. 2. Of particular note is the stronger source practically in the cluster center at 
the position of stars Al-f A2FA3-FB. Its contour lines at A 3 and 6 cm appear slightly elongated 
[major and minor axes of (1 32” ± 0.01”) x (1.09” ± 0.0P ), compared to a beam-size of (0.94 x 
0.77”) at 3 cm, and (1.86” T 0.04”) x (1.28” ± 0.02”) cf. (1.53” x 1.23”) at 6 cm). Its spatially 
integrated flux is S u = 3.3 i 0.3 mJy at 3 cm, as expected if arising in the winds of - 10 - 20 hot 
O/WNL stars due to thermal bremsstrahlung (using the model of Wright & Barlow 1975 together 
with the Kudritzki et al.1987 law, typical wind velocities of 1000 km/s and stellar radii of 10 R©). 
The power law index between 3 and 6 cm is a = -0.2 ± 0.2 (for S„ a z/ a ), significantly smaller 
than a — 0.6 expected for the thermal emission from hot stellar winds. This, along with its 
extended nature, suggests the presence of a non-thermal radio component, possibly arising in zones 
of colliding winds (although without obvious X-ray counterpart). The second central radio source 
is located ~3” West of the lirst source and is only seen at 6 cm, but with nearly the same flux as 
the first source, near stars H49-t-H50. It has dimensions (2.63” ± 0.20”) x (1.52” ± 0.06’) and 
is therefore also extended. Its non-detection at 3 cm strongly suggests non-thermal emission and 
may be caused by colliding winds. Strangely, star C does not appear to be a strong radio source, in 
contrast to its outstanding X-ray flux. This implies that in OB stars the X-ray and radio emission 
mechanisms are not always highly correlated (Waldron et al.1998). 

(7) No statistically significant X-ray variability is seen exceeding 3 a (where a is based on 
Poisson statistics in a given bin) among any of the ten individual, resolved, bright (down to ~100 
counts total) X-ray sources in Fig. la on time scales of 0.5 - 50 ksec. This was determined by 
comparing the standard de\ iation of the X-ray counts about the mean, in 8 equal bins, with their 
instrumental errors, dominated by Poisson statistics. Counts were summed in circles of diameter 
7” outside and 3” inside the central core of Fig. 2. 


3.3. Diffuse X-ray Emission 


Preliminary reduction (adaptive filtering and deconvolution) of the original X-ray image reveals 
the presence of weak extended X-ray emission in the immediate surroundings of the cluster core. To 
enhance the visibility and to quantify this diffuse component, we applied two different approaches. 
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Both approaches attempted to subtract all recognizable point sources using image statistics derived 
in subrasters, with special attention paid to spatial variability of the background level. 

In the first method, each source found by the CIAO wavedetect routine 15 was fit with an ellip- 
tical gaussian model convolved with the Chandra PSF (generated assuming a single-temperature, 
absorbed thermal spectrum with Nh = 1.3 x 10 22 cra" 2 and kT = 0.85 keV) 16 . A fit was per- 
formed using a 32x32 pixel subraster surrounding each source, including a background model of 
the general form: B(x,y) = c+ax-j-by, with a, b and c fitted parameters. In so doing, we assumed 
that large-scale (» 16") background variations prevail, while significant positive fluctuations on 
small scales would have been automatically detected as sources. We held fixed the fit parameters 
of serendipitous sources in the subraster except for normalization, thus preventing these sources 
from biasing the background model. Pixels in the original image were then replaced by the Pois- 
sonian deviations computed from the background model whenever the input image pixel-intensity 
exceeded the background model by 50%. 

In the second approach, we divided the central 510 x 510 pixel part of the image into 9 equal 
subsections, in order to address the issue of gradual variations in the background. We used a 
PSF constructed from a combination of bright, isolated point sources in the original image, with 
stars MTT 68 and MTT 71 as the main constituents. (Pile-up for MTT68 and MTT71 will have 
negligeable impact on the results.) The specially designed source- recognition routine as described 
in §3.2 was then run on the whole 510 x 510 region. 

By slightly varying the detection threshold around the 3 — a level (~10 counts per source) 
we obtained cr-dependent lists of sources, subsequently compiling them via cross-referencing into 
a final list. Then, a properly scaled PSF was subtracted from each point source contained in the 
final list until, after several iterations, the locally calculated (10 x 10 pixel subraster) flux density 
matched the corresponding average value in each of 9 subsections on the original X-ray image. This 
average value was also recalculated for each iteration to acommodate the changes produced by the 
source subtraction. 

Since both approaches yield comparable results, we combined them by simple averaging into 
an image with pixels of size 12" x 12", as shown in Fig. 5. Clearly, the high density of sources in 
the central core reduces the reliability of both techniques there (i.e. inside the region limited by 
a red circle in Fig. 5). Detailed examination of the diffuse emission in the core will be given in a 
forthcoming paper. Therefore, we take advantage of the ~Gaussian form of the spatial extension 
of the diffuse X-ray emission outside the dense core to extrapolate into the core. We thus find that 
the total integrated diffuse X-ray emission, ~ 10 343 erg/s (net unabsorbed, integrated within the 
central, r < 2', region), constitutes ~20% of the integrated point-source stellar X-ray flux. Both 


15 Wavelets are well suited here for identifying structures, including stars, but they are not reliable for obtaining 
stellar X-ray fluxes, as in §3.2. 

16 cf. http://xassist.phys.cmu.edu/ximgfit/manual/nodel6.html 
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this flux ratio, as well as the observed large extension of the diffuse X-ray emission in NGC 3603 
contradict the model of Cani 'pRaga k Rodrigues (2000) applied to similar starburst regions, where 
the extended X-ray emission is predicted to arise from multiple wind- wind interactions. 

To address the possibility that the diffuse emission might be generated by faint point sources 
falling below the 3 -cr detectability level, we obtain a N(S) (number versus source-flux) distribution 
using all the detectable sources, then extrapolate N(S) to lower fluxes. We ignore the fact that the 
form of the distribution is not spatially invariant; N(S) tends to change toward the centre of the 
image. This limits the accuracy of our extrapolation to within an order of magnitude. Otherwise, 
subdivision of the initial image into smaller areas would prevent one from obtaining sufficient 
statistics for constructing a robust N-S diagram for each region. Uncertainty of the point-source 
removal in the very central, overcrowded part of the image (defined approximately by the central 
red circle in Fig. 5) forces us to exclude this region in deriving the N(S) relation. We estimate the 
number of potential point sources in the area occupied by the diffuse X-ray emission as 1.5 x 10 3 , 
which is close to the number of active PMS stars expected in the region. However, such young 
stars may contribute only ^ 10 32 5 erg s" 1 , orders of magnitude below the diffuse background flux. 
If not directly related to PMS stars, the faint undetected point sources would have accounted for 
~ 25% (or even all, allowing for the large uncertainty in the number of extrapolated faint sources) 
of the diffuse background. 

The intensity of the observed diffuse X-ray component does not contradict the non-detection 
of any extended radio emission coming from the same volume of space. Indeed, assuming optically 
thin thermal bremstrahlung emission with kT ~2 keV, a radio flux of ~ io(“^)-(- 5 ) mJy/beam 
at A— 3 cm is expected, which is far below the ATCA 1 o detection limit (0.15 mJy/beam) of our 
observations. 

We have extracted (Fig. 6) the integrated diffuse X-ray component spectrum in an annulus 
around the core region (excluding obvious stellar sources) and a region beyond the core to the 
SE of star MTT68. We fit each region with a single-temperature absorbed thermal model (Mewe, 
Gronenschild k van den Oord 1985; Mewe, Lemen k van den Oord 1986; and Kaastra 1992, 
including the updated Fe-L emission calculations by Liedahl, Osterheld k Goldstein 1995). With 
reduced x 2 = 1*28, the X-ray temperature of the core is kT ~ 3.1 ± 0.3 keV compared to kT « 
2.1 ±0.7 keV in the outer region. The core also had a lower X-ray absorbing column, Nfj ~ 0.7 x 10 22 
cm -2 vs. Nji ^ 2.4 x 10 22 cm -2 in the outer region of the field; this is compatible with the increase 
in E(B — V) going out from the core (Melnick et al. 1989). 

As an additional, independent check on the reality of the weak, extended emission component, 
we performed a census of all the point-like sources by varying the detection limit (namely, 3, 5, 
10- a in accordance with Gehrels 1986) and plot in Fig. 7 the integrated, normalized source flux 
and source number as a function of distance from the tentative centre of the cluster (placed at star 
C for convenience). There is an apparent excess of weak (3 - 5 a) sources toward the centre of the 
cluster at r<l' - 1.5', which can be used as an additional argument backing the reality of the X-ray 
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halo. 


3.4. Star Formation and Proplyd-like Objects 

As found in previous lower- resolution work (de Pree et al.1999), most of the cm radio emission 
comes from extended regions well outside the central core, i.e. mostly anti-correlated with the X- 
ray emission. In particular, mostly thermal radio emission is seen from large numbers of clumps of 
various sizes towards the WSW and SSE from the central core, at projected distances of ~ 30 — 90" 
(1-3 pc) (Fig. la). These are associated with regions of active star formation, with the most 
intense near-IR/radio flux coming from two giant pillar-head ionization fronts at ~ 35" from the 
cluster core. These pillar heads are also sources of OH and H20-masers. Inside this radius, there 
is very little gaseous emission of any kind, except for a diffuse X-ray component from near the core 
where the hot stellar winds are emerging (Fig. 5). The size of this wind-blown cavity suggests an 
age inferior to 1 Myr (Drissen 1999), in agreement with recent estimates of the nuclear age of the 
central cluster, e.g. <1.0 Myr from evolution of the massive stars (Pandey et al. 2000), or 0. 3-1.0 
Myr from PMS evolution (Brandi et al. 1999). Near the southern end of the southern giant pillar 
is located an embedded, very young star cluster (E. Grebel, private communication). It appears to 
be a compact IR/radio source that lies between two adjacent Chandra point sources. This will be 
discussed in a future paper. 

In addition to these extended regions, the three proplyd-like structures detected at near-IR 
wavelengths (Brandner et al.2000) are resolved as very strong, extended cm radio sources (see 
Fig.l), with both thermal and non-thermal components. The NGC 3603 proplyd-like objects are 
the most massive, luminous and distant of their kind discovered so far, and the only ones where 
magnetic processes are inferred due to their non-thermal components (Miicke et al.2001). In the 
meantime, a 4th proplyd-like object has been identified in NGC 3603 (Miicke et al. 2002). However, 
it appears to have more the appearance of a mini-pillar. 


4. Integrated Properties and Comparison with Other Starbursts 

The net instrumental background-corrected count rate from the entire ACIS-I observation is 
~ 1 count s" 1 . Assuming a 3 keV thermal spectrum with Nh = 7 x 10 21 cm -2 and a distance of 
7 kpc, the total luminosity in the 0.5-10 keV band of NGC 3603 in the ACIS-I field of view is Lx 
(0.5-10 keV) = 3 x 10 34 erg/s (net observed) or 1 x 10 35 erg/s (net unabsorbed). 

With time the massive stellar population will evolve and eventually explode as supernovae, 
producing X-ray bright supernova remnants and compact objects, so that the integrated X-ray 
luminosity from the region will increase (Van Bever & Vanbeveren 2000). Other older, well studied 
starbursts like 30 Dor in the LMC and M82 show X-ray emission dominated by SNRs and collapsed 
objects. The integrated luminosity of 30 Dor measured by ROSAT is ^ 6 x 10 37 ergs s -1 in the 
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0. 1-2.4 keV band, with only 1-2% of the emission produced by unevolved massive stars (Norci & 
Ogelman 1995). The X-ray temperature of 30 Dor from the ROSAT observations, however, is 
only 0.3 — 0.4 keV, significantly lower than the temperature of NGC 3603 measured by Chandra. 
ASCA observations of the canonical starbursts M82 (Tsuru et al. 1997) and NGC 253 (Ptak et 
al. 1997) showed that these systems contain distinct hard and soft emission components, with 
the hard component being more centrally concentrated. These results were confirmed by Chandra 
(Griffiths et al. 2000) and XA/M-Newton (Pietsch et al. 2001) observations, which found that most 
of the hard component can be; attributed to multiple (often variable) point sources and an extended 
nuclear component (at temperatures on the order of 5-10 keV), whereas the soft component is due 
to disk emission and a superwind outflow. As has been known from ROSAT PSPC and HR I 
observations (Strickland, Ponman & Stevens 1997), the disk/superwind emission in these systems 
tends to have a temperature of the order of ~ 1 keV, decreasing to ~ 0.4 keV in the halo, and 
X-ray enhancements in the vicinity of known HII regions are often observed (e.g., M51 in Marston 
et al. 1995). It is therefore possible that some fraction of the ( Lx ~ 10 40 ergs s *) disk emission 
in starburst galaxies is due tc> < 10 5 NGC 3603- type HII regions, while the bulk of the hot ISM 
outside their nuclei has expanded and cooled to lower temperatures than those found in NGC 3603. 


5. Summary and Conclusions 

The most important findings of this work are: 

(1) While its more luminous member stars reveal a general trend between Lx and LboI as seen 
in other less exotic clusters, the scatter is enormous and definitely intrinsic to the stars. Some of 
the X-ray overluminous sources may indicate the presence of colliding winds in binaries. 

(2) The three WR stars appear anomalously bright in Lx compared to other WR stars. Prob- 
ably, these are not ordinary WR stars, rather, extremely luminous main-sequence stars with winds 
of relatively low X-ray opacity. Another factor that could enhance their Lx flux is colliding-wind 
emission in the two WR stars that are suspected to be binary. 

(3) Large numbers of X- ray sources are found with greater frequency towards the cluster center 
and no obvious optical counterparts. If these sources are cluster members, then they may be active, 
low-mass PMS stars, though other types of high Lx / Lq 0 i objects cannot be entirely ruled out (for 
example some of them might be isolated neutron stars, or heavily embedded X-ray binaries). 

(4) Spatially integrated X-ray spectra reveal that Nh increases along with E(B — V ), as one 
goes from the cluster core to the periphery, as expected for a wind-blown cavity of cooler gas and 
dust evacuated by an ensemble of very hot stars. 

(5) No significant variability of individual X-ray sources was found on time scales ranging from 
0.5 to 50 ksec. 

(6) Central diffuse X-ray emission may arise from multiple merging/colliding hot stellar winds, 
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with negligible contribution from unresolved PMS stars. 

(7) Most of the cm radio emission arises from peripheral regions of current star formation, 
with little or no correlation with X-rays. However, an extended, possibly non-thermal component 
(due perhaps to colliding winds) is also present in the central cluster region. 

(8) NGC 3603 appears to be a somewhat younger and hotter, scaled-down version of typical 
starbursts found in other galaxies. 

In a future paper, we will discuss our recent Chandra images of the SMC’s youngest, dense 
cluster NGC 346 at 0.1 solar abundance (e.g. Massey et al.1989) for comparison with NGC 3603 
at solar abundance. 
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Fig. i _ Overview 3.6' x M.6' (7.3 x 7.3 pc at the adopted distance for NGC 3603 of 7.0 kpc) 
maps in linear-intensity cerilered on NGC 3603 in (a) smoothed X-ray flux (corrected for exposure 
variations and cleaned of detector artifacts), and (b) groundbased optical emission. Both have the 
same scale and resolution ('-■ 1"). Radio contours (blue, resolution also ~ 1", with contour levels 
of 3, 5, 7, 10 -cr, where 1 a - 1.5 x 10~ 4 Jy beam' 1 ) are superposed on Fig. la and b, with the 3 
proplyd-like objects indicated. Identifications (cf. Melnick et al.1989) of member OB stars outside 
the central core and with detected Chandra X-ray fluxes are indicated in Fig. lb (two stars are 
repeated in Fig. la for orientation purposes). The central red boxes indicate the enlarged region 
in Fig. 2. 
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Fig. 2. — Enlarged 21" x 21" (0.71 x 0.71 pc) region centered on NGC 3603, with superposed 
Chandra X-ray (green) and HST optical (red) images. The HST image has been Gaussian- smoothed 
to a resolution of ~ 0.5", to match the deconvolved and Gaussian-smoothed X-ray image. Radio 
contours in steps as in Fig. 1 are superposed. Overlapping (within a PSF radius) X-ray and optical 
sources become bright yellow; these are identified with prefix H, as members of NGC 3603 (Moffat 
et al.1994), with star H7b added. Note that the bright central stars Al, A2, A3, B and C are H30, 
H31, H26, H23, H18 of Moffat et al. (1994). When more than one star coincides within a PSF 
radius, appropriate identification is given. Epoch 2000 equatorial coordinates are indicated. 


Fig. 3. — Composite summed X-ray spectrum of 7 fairly well resolved early-0 stars that do not 
suffer pile-up in the core of NGC 3603, with model superposed. See text for details. 


Fig. 4. — X-ray versus bolometric luminosity of all X-ray detected member hot-stars of NGC 3603 
(labelled stars in Fig. lb: circles and Fig. 2: triangles). MK spectral types are indicated when 
available, from Drissen et al. (1995) and Skalkowski et al. (in prep.). The dash-dotted line follows 
the relation Lx/ L boI — 10 -7 . The horizontal dashed line shows the 3-cr detection limit for stars on 
the ACIS-I image. We make no attempt to plot individual stars below this limit. Since no correction 
for pile-up has been applied, the derived X-ray luminosities for stars brighter than Lx 10 34 erg 
s" 1 are underestimated roughly as indicated by the arrows. Error bars are not visible when they 
are smaller than the size of the symbol. 


Fig. 5.— Grayscale: rebinned and smoothed central part of the Chandra/ACIS-I image. All point 
sources (in green, from the raw X-ray image) were removed to reveal the extended, diffuse X-ray 
component. Two point sources are labeled in accordance with Fig. 1. 


Fig. 6. — The integrated diffuse component of X-ray surface brightness spectra at two locations 
near NGC 3603. The upper plot is the spectrum in a circular annulus centered near star C at 
the center of the cluster and extending from 7" to 22", excluding obvious stellar sources. The 
lower plot shows the spectrum in a circular region of radius 12.3" centered to the SE of star 68 at 
a = 11 : 15 : 02, 6 = —61 : 14 : 53. The solid histograms are the best fit models as described in the 
text. 


Fig. 7. — The integrated, normalized source flux (upper panel) and source number (lower panel) as 
a function of the distance from the approximate centre of the cluster (star C). Full line shows all 
the sources detected at > 3 a level; dotted line corresponds to > 5cr; dashed line shows the > 10a 


sources. 
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